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ABSTRACT 

(  In  this  paper,  the> expressions  of  plume  axis  height,  reflecting  coefficient 
and  transporting  factor  in  Gaussian  Diffusion-Deposition  Mode!  are  derived 
by  using  theoretical  and  experimental  results  in  planetary  boundary  layer  in 
order  to  obtain  a  new  Gaussian  Diffusion-Deposition  Model  that  could  be 
used  to  describe  variation  of  wind  speed  with  height. 

The  effects  of  wind  speed  profile  on  relative  ground-level  concentrations 
are  also  discussed  by  using  the  calculated  results  obtained. 

I .  FOREWORD 

The  diffusion  theory  based  on  the  Gaussian  model  and  the 
results  of  the  experimental  studies  have  received  wide  attention 
among  the  scientific  community.  In  the  study  of  diffusion  of 
particulate  matter,  the  Gaussian  diffusion-deposition  model  that 
evolved  from  considerations  of  the  gravitational  settling  of  mass 
points  and  the  effect  of  reflection  from  the  ground  surface  has 
been  widely  applied  because  it  still  retains  the  fundamental  proper¬ 
ties  of  the  Gaussian  model  [1].  In  this  model,  the  ground  level 
concentration  assumes  the  expression 

C(*,  y,  0,  H)  -  exp  { -  (y72*}  +  (Jf-  7^/t.) 72*3  > 

a’(J)*1  ~[l^TV;+(^JT-k)^(dcydx)  J 

In  the  above  equations,  Q  is  the  intensity  of  the  source,  H  is  the 

height  of  the  source,  C  is  concentration,  and  a z  are  diffusion 

parameters,  u  is  wind  speed,  V  is  gravitational  settling  velocity, 

s 

is  deposition  velocity,  x,  y  and  z  are  coordinates  and  otQ(x)  is 
the  reflection  coefficient  of  the  ground. 

It  is  worth  pointing  out  that  in  the  Gaussian  model,  the  average 
wind  speed  acts  as  a  constant,  and  has  a  transporting  effect  on  the 


diffusing  mass  points.  Here  the  relation  between  the  concentration 
C  and  the  wind  speed  u  is  given  by  Cctu-1.  On  the  other  hand,  in 
the  Gaussian  diffusion-deposition  model,  the  effect  of  the  wind 
speed  is  threefold.  Not  only  does  the  wind  speed  have  a  transport¬ 
ing  effect  on  the  pollutants,  but  it  also  has  an  effect  on  the  traj¬ 
ectory  of  mean  motion  of  the  diffusing  mass  points  and  the  reflect¬ 
ion  coefficient.  Moreover,  since  in  fact  the  wind  speed  varies 
continually  with  height,  the  assumption  of  a  constant  wind  speed 
necessarily  produces  some  error  in  the  estimation  of  concentrations 
in  the  downwind  direction. 

Wojciechowski  [2]  discussed  the  effects  on  the  diffusing  par¬ 
ticles  of  the  wind  given  according  to  the  expression  by  Roberts.  In 
this  paper,  the  theory  and  experimental  results  of  wind  speed  profile 
in  planetary  boundary  layer  studies  will  be  used  to  discuss  the 
effects  of  the  average  wind  speed  profile  on  plume  axis  height, 
reflection  coefficient  and  the  transport  of  diffusing  mass  points. 

The  deviations  of  the  estimation  of  concentrations  arising  from 
these  effects  will  also  be  examined. 

II.  AVERAGE  TRAJECTORY  OF  MOTION  OF  DIFFUSING  PARTICLES 

AND  THE  WIND  SPEED  PROFILE  IN  ATMOSPHERIC  BOUNDARY  LAYER  /157 

In  equation  (1^,  the  height  of  the  plume  from  the  ground  is 
described  with  H-V  x/u.  If  the  wind  speed  is  a  function  of  height, 
then  this  has  an-~effect  on  the  diffusing  particles  with  settling 
velocity  Vg.  The  tra jectr>s^of^such  a  particle  can  be  described  by 
the  following  equations: 

/  dx 

f  • 

_ii~  (2) 

dt  • 

In  the  equations,  t  is  time  and  u(z)  is  the  wind  speed  at  height  z. 

Given  a  definite  functional  form  for  u(z)  and  suitable  boundary 
conditions,  one  can  obtain  the  equation  for  the  trajectory  of  the 
motion  of  the  particle. 
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From  the  similarity  hypothesis  applied  to  the  ground  level 
boundary  layer,  the  following  definite  expression  has  been  derived 
for  the  wind  speed  profile  [3,1]: 


Under 

neutral  conditions, 

• 

stable  conditions, 

«(*)  -(«•/*)  In(*/*t) 

(3) 

Under 

Under 

unstable  conditions. 

n (»)  -  (u./i)^ ln-1-  +  ?-£-) 

(4) 

«(.)  -  («./*)  (in-: - 2  ]„  -  In  _ 

l  *.  2  2 


(5) 


+  2tan-'(i-i8*/7,)l/‘  — j-J 

In  the  above  equations,  u*  is  frictional  velocity,  z^  is  roughness 
length,  L  is  Monin-Obukhov  length,  k  is  von  Karman's  constant,  and 
3  is  a  constant  to  be  determined. 


In  fact,  as  the  thickness  of  the  constant  stress  layer  is 
usually  around  several  tens  of  meters,  the  above  results  have 
limited  applications.  In  diffusion  studies,  the  power  law  express¬ 
ion  is  often  used  instead,  i.e., 

«(*)  —»,(*/*,)'  (6) 

In  the  above  equation,  u^  and  are  reference  velocity  and  height, 
and  the  index  P  is  a  parameter  that  varies  with  degree  of  stability 
and  degree  of  roughness  of  the  ground. 

If  the  initial  conditions  are 

<«o,  *-«, 

the  wind  speed  profile  is  taken  separately  as  that  given  in  equations 
(3)-(6),  equation  (2)  is  solved  and  let  z  (plume  axis  height).  =  z, 
then  one  can  obtain  the  equations  for  the  trajectory  of  the  diffus¬ 
ing  particle  as  follows: 

(7) 
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v*« 24^(2  In  1/  —  2  —  2  In  ». +  6  !»  2  -*) +  3/)(l  -  16-tf/£)  (2  !n(l  -  ibH/Ly 

-  ln(  1  —  1  tfltLy*  -  1  -  2  tan“‘(l  -  llH/Ly*  )  +  67' tan"'  0  “  1  •')'?* 

+  2t(l-16 U/T.)"4  -6^(1 -16 H/L)ut 

C-20(2lnT-2-2ln«,  +  Sln2-^  +  3la-16  V/OC2!n(l-16_»A)l/4 

-  ln(l  - 16  »/£)*'*  -1-2  tan-‘(-l  - 16  »/£)U4 )  +  6-&  tan“*Cl  - 16  »/■£)*'* 

1  +27,(1- 16  */^)W4  —  6£Cl  —  15"»/-t)1/<  (9) 

I -I7fl--^(i,+  l)l7rn  (10) 

- - fe 

Equations  (7) -(9)  give  the  relation  between  plume  axis  height  and 
distance  under  neutral,  stable  and  unstable  conditions  in  an 
implicit  form,  while  equation  (10)  gives  in  an  explicit  form  the 
relation  between  the  plume  axis  and  the  downwind  distance  for  the 
case  where  the  wind  speed  has  an  exponential  distribution. 

On  the  basis  of  Pasquill  and  Smith's  work  [4],  for  flat  lands 
of  the  open  countryside,  the  relation  between  the  Monin-Obukhov 
length  and  Pasquill' s  stability  classification  is  as  shown  in  Table 
1.  The  values  for  P  for  different  degrees  of  ground  roughness  and 
stability  are  taken  to  be  those  given  in  Table  2.  Thus,  based  on 
various  roughness  lengths,  one  can  calculate  for  fixed  values  of 
Vg/u  or  the  height  of  the  plume  axis  z  under  various  stability 

conditions . 

Consulting  Table  2,  we  chose  two  sets  of  values  for  P  to 
represent,  respectively,  the  various  stability  conditions  of  cities 
and  countryside.  We  also  took  u>  =  0.5,  0.1,  0.01  and  H  =  100m  to 
perform  calculations  on  equation  (10) .  The  results  are  shown  in 
Figure  1.  It  can  be  clearly  seen  from  this  figure  that  when  P  =  0, 
indicating  that  the  wind  speed  is  constant,  equation  (10)  transforms 
into  the  form  given  in  equation  (1) ,  and  the  trajectory  of  motion 
of  the  particle  becomes  a  straight  line  intersecting  the  x-axis  at 
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T*bl«  1  Ralitioaa  between  Puquill  clM«fie»tioa  ond  Monin-Obaibor  length  L 
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1 — roughness  condition;  2--flat  countryside;  3 — cities  (Japan) ; 

4 — western  suburb  of  Peiching;  5 — Shenyang;  6 — references;  7 — remark; 
8 — Remark;  This  is  taken  from  an  evaluation  of  environmental 
quality  of  the  western  suburb  of  Peiching 


a  slant.  P  /  0  indicates  that  the  wind  speed  varies  with  height, 
causing  the  plume  axis  to  bend  downwards.  With  increasing  stability 
(increasing  value  of  P) ,  the  curvature  of  the  curve  increases  and 
the  deviation  of  the  trajectory  from  a  straight  line  becomes  more 
pronounced . 

Applying  the  same  method  to  perform  calculations  on  equations 
(7) -(9)  for  given  values  of  Vg/u*,  zQ,  L,  3  and  other  related  para¬ 
meters,  one  can  obtain  results  similar  to  that  given  in  equation 
(10)  . 

Figure  2  gives  the  results  obtained  from  equations  (7)  and  (10)  , 
under  neutral  conditions  and  for  roughness  lengths  z^  =  1  cm  and  50 
cm.  It  shows  that  the  results  obtained  from  these  two  equations  are 
very  similar  wlieh  the  roughness  length  is  small. 
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Pig.  2.  Comparison  between  trajectories  of  mean  motion  of 
particles  in  neutral  conditions  when  wind  speed  profile 
is  shown  respectively  by  power  law  or  logarithmic  law. 


1 — neutral  conditions;  2 — power  law;  3 — logarithmic  law 


III.  GROUND  REFLECTION  COEFFICIENT,  TRANSPORTING  FACTOR 
AND  WIND  SPEED  PROFILE 


The  ground  reflection  coefficient  in  equation  (1),  a  (x)  ,  has 

been  derived  under  the  assumption  that  the  wind  speed  u=uu  [7] . 

n 

If  the  variation  of  wind  speed  with  height  takes  the  form  given 
in  equation  (6) ,  then  the  plume  axis  will  be  as  shown  by  the 
broken  lines  in  Figure  3.  If  the  streamline  for  constant  wind  speed 
is  represented  by  ST,  then,  because  the  wind  speed  increases  with 
height,  the  streamline  moves  backward  to  ST'.  The  ground  reflect¬ 
ion  coefficient  a^(x)  can  be  derived  as  follows. 

The  streamline  ST'  should  be  described  by  the  equal  i  below: 


s 


Fig.  3.  Sketch  of  effects  of  wind  speed 
with  height  on  plane  axis. 


I  -o>ui+i)zr/n)^T  (id 

*.(*!■) 

„  ★ 

In  the  above  equations,  z  denotes  the  distance  from  ground  level 
of  the  particle  in  the  plume  on  the  streamline  ST' .  Differentiat¬ 
ing  equation  (11),  we  obtain 


Cl— <P+l)*r 


(12) 
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Equation  (12)  shows  that  the  vertical  velocity  of  the  particle  is 
made  up  of  the  gravitational  settling  velocity  and  the  turbulent 
flow  diffusion  velocity.  The  negative  sign  denotes  downward  trans¬ 
port  . 


Define  the  rate  of  deposition  D  as  the  product  of  the  deposi¬ 
tion  velocity  and  the  ground  level  concentration 

D  =  VdC(z,y,  0)  (13) 


If  deposition  occurs  at  ground  level  and  ct^(x)  is  used  to  denote 
the  consumption  of  the  source,  then,  according  to  the  law  of  con¬ 
servation  of  mass,  the  ground  level  settling  should  be  equal  to  the 
difference  between  the  fluxes  from  the  real  and  image  sources. 
Hence, 


V<(1  +«,)<?.(«,  y.  Q)  -(F.+  V,)C,(x,  y,  0)  -  a,  (V«  +  V,)C, (*,  y,  0) 


(14) 


In  the  above  equation,  Cg  denotes  that  part  of  the  ground  level  con¬ 
centration  due  to  the  real  source;  Vfc  denotes  the  turbulent  flow 
diffusion  velocity.  From  equation  (12),  on  the  ground  level 

Inserting  equations  (12)  into  (14), 
one  obtains  the  expression  for  a^(x): 


-{ 


2  Vt 


r4+v,+ J 


(15) 


In  the  Gaussian  diffusion-deposition  model,  the  wind  speed 
used  in  describing  the  transport  effect  of  the  wind  is  usually 
taken  to  be  that  at  the  height  of  the  exhaust,  u„,  in  the  same 
manner  as  when  using  the  Gaussian  model.  Since  the  height  of  the 
plume  axis  continually  decreases,  taking  the  wind  speed  at  the  height 
of  the  plume  axis  is  a  more  nearly  true  representation  of  the  actual 
situation.  If  the  distribution  of  the  wind  speed  is  given  by 
equation  (6) ,  then  the  wind  speed  at  the  height  of  the  plume  axis 
should  be 


u(  t  ) 


l-P+l 


(16) 
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IV.  EFFECTS  ON  GROUND  LEVEL  CONCENTRATION 


To  help  illustrate  the  problem,  we  use  equation  (1)  as  our 


basis  and  take  u  =  u 


H* 


The  relative  ground  level  concentration 


can  be  expressed  as 

0.0)_(i+qt(ic)) 

2^,0. 


(-(#- Vjs/*Kyfio*) 


(17) 


The  effects  of  the  wind  speed  profile  on  the  ground  level  concentra¬ 
tion  as  produced  via  plume  axis  height,  reflection  coefficient  and 
the  transport  term  can  be  obtained  by  calculating  C. ,  C2  and  C^  by 
substituting  equations  (10),  (15)  and  (16)  into  the  corresponding 
terms  in  equation  (17).  Let  I ^ ,  I  and  1^  denote  the  relative  change 
in  concentration  due  to  the  above  three  factors,  and  we  have 


r  n,r>  r  I*- (ff -<**)*') 

J.-C/C.-exp^ - ^ - L  J 

r.-c/c.-fi-: - _ V 

L  V*  +  V.  +  (  u,g  -  v.x) 


(18) 


(19) 


,  P 

/,-C/C,-[i  -  <1»+ 1) 


(20) 


After  the  variation  of  wind  speed  with  height  is  taken  into  consider¬ 
ation,  the  axial  relative  ground  level  concentration  takes  the 
expression 


«xC(x,  o,  0) 
- ? - 


_i±f«^L.exp(-W20 


(21) 


For  a)  =  0.1,  V,  =  V  =  0.1  m/sec,  H  =  100m,  and  a  and  a  taken 

CL  S  y  2 

from  the  interpolation  formula  of  Briggs  [6],  perform  calculations 
on  equations  (17) -(21)  for  a  given  value  of  P .  As  the  value  of 
u(z)  in  equation  (21)  self-adjusts  with  plume  axis  height,  we  set 
the  lower  limit  of  the  height  to  be  u(z)  =  u^g . 
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Pig.  4.  Variations  of  concentration  along  relative 
plume  axis  with  distances. 
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Table  3  gives  the  variation  of  1^,  I ^  and  with  distance. 

It  can  be  seen  from  the  table  that  1^  first  increases  then  decreases 

with  distance.  In  the  range  of  1-4  km,  1^  >  1,  and  I^  and  I^ 

increase  with  increasing  distance.  However,  is  greater  than  1 

while  I ^  is  less  than  1.  The  rate  of  variation  of  I, ,  I^  and  I^ 

with  distance  increases  with  increasing  stability.  This  shows  that 

the  assurrmtion  of  u  =  u„  =  constant,  through  its  effect  on  the 

n 

plume  axis,  has  caused  the  estimation  of  the  concentration  in  the 
range  of  1-4  km  to  be  on  the  high  side.  The  same  assumption, 
through  its  effect  on  the  ground  reflection  coefficient  and  the 
transport  factor  has,  respectively,  led  to  an  overestimation  and 
underestimation  of  the  axial  ground  level  relative  concentration, 
the  effect  of  the  wind  speed  profile  on  concentration  being  greater 
for  higher  stabilities. 

The  results  computed  from  equations  (17)  and  (21)  are  shown 
in  Figure  4.  It  can  be  seen  that  for  a  roughness  length  of  1  cm, 
and  type  C  atmospheric  stability,  the  results  obtained  from  equa¬ 
tions  (17)  and  (21)  are  relatively  close.  Out  beyond  3-4  km,  the 
results  obtained  from  equation  (17)  are  lower  than  those  obtained 
from  equation  (22) .  Under  neutral  conditions  (type  D) ,  the  differ¬ 
ence  becomes  more  pronounced.  Under  stable  conditions  (types  E  and 
F) ,  the  ground  level  relative  axial  concentration  calculated  from 
aquation  (21)  has  an  increased  rate  of  change  with  distance.  This 
shows  that  the  overall  effect  of  the  variation  of  wind  speed  with 
height  on  ground  level  concentration  is  an  increased  maximum  value 
of  concentration  and  the  appearance  of  the  phenomenon  of  a  "regress¬ 
ion"  of  distance.  This  result  is  in  keeping  with  the  observations 
made  by  Csanady  in  1964  [9].  The  fact  that  the  difference  between 
the  results  obtained  from  equations  (17)  and  (21)  increases  with 
increased  stability  indicates  that  the  effect  of  wind  speed  on 
ground  level  relative  concentration  is  more  pronounced  under  stable 
conditions,  and  should  be  paid  attention  to. 
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V.  CONCLUSION 


From  the  above  analysis,  it  can  be  seen  that  the  newly  derived 
Gaussian  diffusion-deposition  model  that  takes  into  account  the 
variation  of  wind  speed  with  height  is  a  better  description  of  the 
diffusion  of  particles  in  the  atmosphere.  The  wind  speed  profile 
can  affect  the  ground  level  concentration  through  curving  of  the 
plume  axis,  the  reflection  coefficient  and  the  transporting  factor, 
and  the  effect  becomes  more  pronounced  under  the  stable  layer. 

(Received  on  January  8,  1981) 
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